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DECUS Program Library Write-up 


DECUS No. 8-283 


Summary 

This program computes the means, between-groups sum of squares and mean square, and within- 
groups sum of squares, mean square, and standard error for univariate data. Versions of the 

program are available for carrying out transformations on the data on input i e • 


• 

angular transformation 

(arcs in N/~x ) 


square root transformation 

(\/ x + 0.5) 


logarithmic transformation 

t lo 9| 0 (x + 1) _7 


Tapes Required 


Form of p rogram tape - The program is written in the PDP-8 FORTRAN-D 
source language. 


language, and is in the 


Y l* 711 i°f ^ q / a analysed should be punched on to paper tape in 

T e data should begin with the number of groups, followed by the number of values 
group and the values of this group, the number of values in the second group, e.g. 


the ASCII code, 
in the first 


4 

5 

40 24 46 20 35 
5 

29 27 20 39 45 
5 

11 31 17 37 39 
5 

17 21 28 33 21 


Operating Instructions 


• FORT 

*OUT-S:AVSC 

* 

*IN-R: 

* ^ 

t 

•READY 

T 
















If the program has already been compiled on to the disk the program may be called back into 
core as follows: 

. FOSL 

*IN-S:AVSC 

* 

*OPT- 

* ^ 

*READY 

t 

Output 

The program prints the mean for each group, together with the within-groups degrees of freedom 
sum of squares, and mean square, and the standard error of the group means. These are followed 
by the mean of all groups, the between-groups degrees of freedom, sum of squares, mean square ^ 
and F-ratio. The results for the example of Form of data tape are given below in the form in 
which they are printed. 

0.330000E+2 
0.320000E+2 
0.270000E+2 

0.240000E+2 16 0.164800E+4 0.103000E+3 0.453872E+1 

0.290000E+2 3 0.270000E+3 0.900000E+2 0.873786E+0 

Storage and Limitations 

Normal for FORTRAN-D. The maximum number of individuals in a single group has been set at 
100/ but there is no restriction on the number of groups. 


Method 


The method of analysis follows closely the conventional method of calculating the analysis of 
variance for a single classification. 
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AVSC 


sisa* 


L 

C 

1 

100 


20 

101 


60 


10 2 


9 


103 


UARian,( 

va*vuv^ u y ci00 > 

2.l00^ a,V9a « 
FORMAT Cl) 2 

*> 9 *=1 J 8 

v^ 100 '* 3 

20 J=D,M3 
vJjL 2 . J r 101 ' VC J) 

co\itimue F( V(J) + 1 - 0 >* 
FORMAT ( E, £ E „ 
V4=V5=0. 0 
DO 60 M 4= 1 , 3 

U 4 =V 4 +VCM 4 ) 

CMmuE < ' l4> * V< '' 4,> 

V2= V2+ V4 
V3=V3+V5 

V 6 = V 4 / VI 0 


T/PET 102, V6 
format C/,E) 

V1 = V 1 + VI 0 
V 4= V 4* V 4/ V 10 

1 / 0 = V 0 + V 4 

V5=V5-V4 
V9=V9+ V 5 
COMTI 1UE 
M 5= V 1 


M 5= iM 5- ,\J 2 
V6=M5 
V7= V9/V6 
V4= \I 2 

V 6 = SQTF( V7/C V 1 /V 4 ) ) 
TK'PE 103,M5, V9, V7> V 6 
FORMAT (I, E, E, E) 

V3 = V2/V1 
M 2= M 2- 1 

V0= V0- < V2* V2/ VI ) 


V 4 = M 2 


V2= V0/ V4 



fl> Vo i 


\|2, V0»V2>V3 

£,1i£»E , e,/ ’ 


-■» AVSe 
0 


• 43429 4 


/>/>/>/) 


““ ffirisr"”" 


190 


20 

101 


60 


102 


103 


cxs 3 F VAHIA^C 6 -’ A ysC 

°if v 'i?v2--«= , 4=^= a - 0 

READ 2»' ^ ’ 

F 3R« t »» 

n3 9 i - 1 > ^ ^ 

READ 2, 100^ 3 
V/10^3 

o: » 2 ® J / 0 v(J) 

READ 2 

CyMTHUE £) 

F1RMAT CE, E> E> *> E) 
V4=V5=0*^ 

DO 60 M4*l* > * 3 

CD OTI OUE 
V2= V2+V4 
V 3= V3+V5 
V6=V^/V10 
TYPE 102* V6 
F D RMAT C/,E) 

V1 = V1 + V1 9 
V4=V4* V4/ VI 0 

V0=V0+ V4 
V 5= V 5- V 4 
V9 = V9 + V5 

cdoti OUE 
0 5= VI 
0 5= 0 5- 0 2 
V6= si 5 
V7= V9/ V6 
V4= 0 2 

V6= SQTFCV7/C VI/V4> > 

TfPE 1 0 3* 0 5* V9 > V7» V 6 
FORMAT (I*E*E>E) 

Mi= V2/ VI 
0 2= 0 2- 1 

V0= V0- C V2* V2/ VI ) 

V 4=0 2 
V2=V0/V4 
V3= V2/V7 

T/PE 10 4* \/S> 02* V0> V2* V3 

I. E.E.E./, 

bj Tj 1 


* 




